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A SEARCH FOR LOST PLANETS IN THE KEPLER MULTI-PLANET SYSTEMS AND THE 


DISCOVERY OF A LONG PERIOD, NEPTUNE-SIZED EXOPLANET KEPLER-150 F 


Josepu R. Scumitt!, JoN M. JENKINS,”, DEBRA A. FISCHER! 


ABSTRACT 


The vast majority of the 4700 confirmed planets and planet candidates discovered by the 
Kepler space telescope were first found by the Kepler pipeline. In the pipeline, after a transit 
signal is found, all data points associated with those transits are removed, creating a “Swiss 
cheese”-like light curve full of holes, which is then used for subsequent transit searches. 
These holes could render an additional planet undetectable (or “lost”). We examine a 
sample of 114 stars with 3+ confirmed planets to see the effect that this “Swiss cheesing” 
may have. A simulation determined that the probability that a transiting planet is lost due 
to the transit masking is low, but non-neglible, reaching a plateau at ~ 3.3% lost in the 
period range of P = 400 — 500 days. We then model the transits in all quarters of each star 
and subtract out the transit signals, restoring the in-transit data points, and use the Kepler 
pipeline to search the transit-subtracted (i.e., transit-cleaned) light curves. However, the 
pipeline did not discover any credible new transit signals. This demonstrates the validity and 
robustness of the Kepler pipeline’s choice to use transit masking over transit subtraction. 
However, a follow-up visual search through all the transit-subtracted data, which allows for 
easier visual identification of new transits, revealed the existence of a new, Neptune-sized 
exoplanet. Kepler-150 f (P = 637.2 days, Rp = 3.86 Ra) is confirmed using a combination 
of false positive probability analysis, transit duration analysis, and the planet multiplicity 


argument. 
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1. INTRODUCTION 


The Kepler space telescope (Borucki et al. 2010) 
has discovered more than 4700 Kepler Objects of In- 
terest (KOIs) that are classified either as confirmed 
planets (~ 2300, e.g., Rowe et al. 2014; Morton et al. 
2016) or planet candidates (~ 2400, e.g., Cough- 
lin et al. 2016), according to the NASA Exoplanet 
Archive (Akeson et al. 2013). The vast majority 
of these confirmed planets (CPs) and planet can- 
didates (PCs) were initially discovered through the 
Kepler pipeline. One of the first steps of the pipeline 
is the Transit Planet Search (TPS) algorithm (Jenk- 
ins et al. 2010a). This searches the entire Kepler 
data set for potential transit signals, called Thresh- 
old Crossing Events (TCEs, Tenenbaum et al. 2012, 
2013, 2014; Seader et al. 2015). Follow-up vetting 
for potential PC status is then performed, the latest 
version of which is described thoroughly in Section 
3 of Coughlin et al. (2016). Additional vetting, such 
as follow-up observations, detailed light curve anal- 
ysis, or statistical methods, can be used to confirm 
their planetary status or rule them out as false pos- 
itives (FPs). Many of these KOIs are in systems 
with multiple KOIs, which allow for easier statis- 
tical validation of their planetary status (Lissauer 
et al. 2012). As such, about half of the CPs dis- 
covered with Kepler data are located in confirmed 
multiple planet systems. 
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A detail in the TPS algorithm is its treatment of 
targets with multiple detections of TCEs. If a TCE 
is found for a target, TPS removes all the TCE’s in- 
transit data points before TPS is rerun on the same 
light curve to search for more signals, creating what 
the Kepler team call a “Swiss cheese” light curve 
(Twicken et al. 2016). This masking of data points 
can hide the existence of additional planets whose 
transits overlap with the previously found TCEs. 

An additional wrinkle arises when trying to dis- 
cover long period planets in multiple planet sys- 
tems. Large, long period planets with few tran- 
sits are often best found by visually inspecting the 
light curves, especially for planets with < 3 transits, 
even more so if they only transit once (Wang et al. 
2015; Uehara et al. 2016). In fact, a citizen science 
program called Planet Hunters (Fischer et al. 2012) 
that allows people online to visually search the Ke- 
pler data for exoplanets specializes in finding long 
period planets. Through the power of visual in- 
spection, Planet Hunters has led to the discovery of 
three exoplanets (Schwamb et al. 2013; Wang et al. 
2013; Schmitt et al. 2014b) and nearly 100 exoplanet 
candidates (Lintott et al. 2013; Wang et al. 2015; 
Schmitt et al. 2016). However, while these planets 
might be easy to spot visually, the fact that the light 
curve is filled with so many other planet transits can 
make them difficult to identify as new planets. They 
can be easily mistaken for or assumed to be a tran- 
sit from another, known planet in the system. This 
problem could be fixed if the light curves had all 
known transit signals subtracted out, leaving only 
the previously undiscovered transits in the data. 
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In this paper, we examine a large fraction of the 
systems with three or more CPs and perform three 
separate analyses on them. First, we simulate the 
percentage of true planets missed by TPS because 
of the algorithm’s removal of known in-transit data 
points. In our second test, we attempt to extract 
potential new planets in the data by subtracting out 
the transit signals of the known KOlIs rather than 
masking out their transits altogether, after which 
we then rerun TPS on the new, transit-subtracted 
light curves. Lastly, we then examine the transit- 
subtracted light curves visually to search for evi- 
dence of additional planets. 


2. SIMULATING FOR LOST PLANETS 


Flattening light curves, re-fitting for planets, and 
then re-running TPS is a time intensive process. We 
also expect a higher rate of missed planets for sys- 
tems with more known planets, as this corresponds 
to more data points being removed from subsequent 
transit searches. For these two reasons, we limited 
our study to only systems with three or more CPs, 
according to the NASA Exoplanet Archive (Akeson 
et al. 2013) as of January 6, 2016. Some of these 
systems, however, were removed from the analysis 
in later steps (see Section 3). The final sample in- 
cludes 114 stars (see Table 2). 

For each target star containing three or more 
CPs, we downloaded the long cadence data from the 
Barbara A. Mikulski Archive for Space Telescopes 
(MAST). In this sample of stars, the maximum time 
baseline was 1470 days. For each of the 114 stars in 
the sample, we then injected a planet into the light 
curve. The planet’s period and the star’s mass and 
radius were used to calculate the transit duration 
(assuming inclination i = 0), and a random epoch 
was chosen. The simulated planet’s transits were 
counted as detectable if at least 50% of the transit 
was contained within the data (ie., not in a data 
gap). The total number of transits detectable for 
each injection was then recorded. For this purpose, 
only the duration of the transit mattered, not the 
depth. We then repeated this for many periods and 
epochs, injecting 1000 planets into each one day pe- 
riod bin from 2 to 1472 days uniformly distributed 
in period and phase (e.g., 1000 planets with a period 
between 2 and 3 days, 1000 planets with a period 
between 3 and 4 days, etc.). A total of 169,050,000 
planets were simulated. 

This procedure was then repeated (with the same 
simulated planets) on the light curve after all in- 
transit data points from known KOIs were removed. 
Only those in-transit data points that we were able 
to successfully remove in our subsequent analysis 
(see Section 3) were removed in this step. Therefore, 
a small number of KOIs did not have their in-transit 
data points removed at this point. 

The removal of the in-transit data points of the 
KOlIs resulted in some of the transits that were 
originally detectable (pre-removal) becoming unde- 
tectable. This change in the window function (the 
Swiss cheesing of the light curve) generally would 
not be a problem when the planet transits the star 
many times. However, TPS requires three transits 


to register a detection, and the loss of one or more 
transits may bring a planet that had 3+4 transits 
below that threshold. These planets are then no 
longer detectable. These are the “lost planets”, the 
planets that would have been detected if in-transit 
data points were properly corrected for instead of 
removed. 

There exists two shortcomings of this simulation. 
One is that, for planets with five or more tran- 
sits, the removal of 2+ transits in a certain way 
could cause the subsequent planet detection to have 
an alias of the true period. For example, remov- 
ing the second and fourth transit in a system with 
five consecutive transits results in a detected pe- 
riod double that of the true period. However, such 
events are rare, so its effects on the period detection 
are therefore ignored. Another complication is that 
this simulation did not test to see how the transit 
signal-to-noise was affected, only how the number 
of detectable transits was changed. Removing data 
points could reduce the signal-to-noise of undiscov- 
ered transits below the detection threshold. This 
implies that we are underestimating the number of 
lost planets caused by transit masking. See Sec- 
tion 6 for a more detailed discussion. 

The probability of originally detectable planets 
becoming undetectable after the removal of the 
known in-transit data points is shown in Figure 1. 
The black line is the probability averaged over all 
stars, while the transparent gray lines in the back- 
ground are the star-by-star probabilities, so that 
darker areas correspond to higher density regions. 
For the 114 stars in our sample, the lost planets are 
broadly distributed in period in a range of approx- 
imately 200-700 days. (A small number of planets 
are lost below P = 200 days, but this is almost 
exclusively for a small subset of stars that were ob- 
served for fewer quarters than the rest.) The dis- 
tribution plateaus between 400 and 500 days with a 
peak in the 480-500 day period range. This 400-500 
day peak has an average lost planet value of 3.3% 
level, meaning that 3.3% of the observable, tran- 
siting planets in this region would be expected to 
be undetected (“lost”) after removing the in-transit 
points of previously found planets. In other words, 
if the planetary system had a planet in that pe- 
riod range, and if that planet transited, there would 
be a 3.3% chance that it would not have been de- 
tectable due solely to the Swiss cheesing of the light 
curve. The average’s maximum of 4.6% occurs at 
P = 493 days. The star-by-star peaks, on the other 
hand, vary between 3.8-10.0%, with the peaks’ loca- 
tions ranging from P = 135 days to P = 497 days. 
These numbers, while small, are not negligible, and 
therefore imply that there may be a small number of 
missing exoplanets in the Kepler data caused by the 
removal of in-transit data points of known planets. 


3. SEARCHING FOR LOST PLANETS 


We began the analysis with the original Pre- 
search Data Conditioning Simple Aperture Photom- 
etry (PDCSAP, Stumpe et al. 2012; Smith et al. 
2012) fluxes for every Kepler star with 3+ CPs. We 
then removed the stellar variability using the PyKE 
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Figure 1. The probability of a planet that was originally detectable (i.e., had 3+ detectable transits) that became undetectable 
(i.e., had < 3 detectable transits) after the in-transit data points of the successfully fit KOIs were removed. A transit was ruled 
“detectable” if at least 50% of it was contained within the data (i.e., < 50% in a data gap). The probability of losing a planet 
at a certain period averaged over all stars is highlighted in black, while the star-by-star probabilities are shown in a transparent 


gray, so that darker areas correspond to higher density regions. 


kepflatten command in PyRAF (Still & Barclay 
2012). This command divides the light curve into 
chunks (or steps). The step size is usually on the or- 
der of one to a few days. It then fits the light curve 
in a window around (and including) the step with a 
polynomial (ignoring outliers). The window size is 
usually approximately double that of the step size 
so that the edge effects of fitting do not affect the 
center portion of the window. The window and step 
sizes were changed to fit each quarter of each star in 
order to remove as much of the stellar variability as 
possible without removing the transits. These de- 
trended light curves were then stitched together into 
one FITS file. This was often successful as deter- 
mined by eye (see Figure 2). However, a significant 
minority of stars had small, residual variability that 
could not be removed without disrupting the transit 
signal. Attempting to completely fit the stellar vari- 
ability would cause the transits to become partially 
filled in because the in-transit data points would 
not register as outliers. The most egregious cases 
were those in which the frequency and magnitude 
of the stellar variability were greater than or ap- 
proximately equal to the duration and depth of the 
planet, respectively. For these, the fitting procedure 
was unable to adequately fit the stellar variability 
without also removing the transit signal. The worst 
cases were removed from the analysis. 

We then used the PyKE command keptransit 
(Still & Barclay 2012) to fit the KOIs in these sys- 
tems across all observed quarters. For each star, 


this was done iteratively starting with the KOI with 
the largest depth. After this KOI was successfully 
fit, the KOI’s transits were subtracted from the de- 
trended light curve. Starting from this new, transit- 
subtracted and detrended light curve, we then per- 
formed the same procedure for the KOI with the 
next largest depth. This was repeated until all KOIs 
were fit and removed from the light curve. Initial 
parameter guesses for the fits were taken from the 
KOI cumuluative list, accessed 2016 Jan. 8, to get 
the most up-to-date (at the time) fit parameters. 
We chose to re-fit the transits rather than use the 
KOI list’s values as fixed values to correct for any 
differences that our flattening made have induced 
in the transits. A sample of this fitting procedure’s 
results is shown for Kepler-253 in Figures 2 and 3.. 

This fitting procedure was successful for the vast 
majority of cases, but not all. If a majority of a 
star’s KOIs could not be successfully fit, it was re- 
moved from the analysis. Several examples of tran- 
sit timing variations (TTVs) are also present in the 
data (Mazeh et al. 2013; Holczer et al. 2016). The 
keptransit command is not equipped to handle 
TTVs, so these were not properly fit. KOIs with 
TTVs from Mazeh et al. (2013) and Holczer et al. 
(2016) were noted, especially those which were vi- 
sually apparent in the fitting results, in order to 
account for them in the later analysis. In the most 
egregious cases, such as Kepler-90 (Cabrera et al. 
2014; Schmitt et al. 2014a), these TTVs were so 
large as to render the entire fit impossible or use- 
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Figure 2. Thirty day portion of the light curve for Kepler-18 (KIC 8644288). Top: The PDCSAP flux, which was our starting 
point for the analysis. Middle: The PDCSAP flux detrended for variability. Bottom: The detrended flux after removing the 
transit signals from the KOIs in Kepler-18. Note that the y-scaling changes in each panel. 


less. These TTV systems were removed. 

The final count of systems that made it through 
all levels of analysis without being wholly removed is 
114 stars, which host 397 CPs and 14 PCs. Of these, 
eight CPs and two PCs around nine stars were not 
successfully fit (see Table 2). These systems were 
still included in the analysis. 


4. TRANSIT-SUBTRACTED SEARCH RESULTS 


We then searched the detrended,  transit- 
subtracted light curve of all 114 stars with TPS us- 
ing (JON: HOW MANY CPU HOURS OR SOME 
OTHER UNIT] on the NASA Pleiades supercom- 
puter. TPS found 33 new, unique signals in 24 stars 
that did not correspond to known KOIs. The new 
signals had between three and six transits. Each of 
these signals were cross-checked with the locations 
of the removed transits of known planets. There 
were 13 new signals that overlapped with at least 
one transit of a KOI that had been subtracted out. 


Two of them overlapped with two removed KOI 
transits, and four overlapped with three removed 
KOI transits. The other seven signals only over- 
lapped with one removed KOI transit. The other 
20 signals did not overlap at all with the removed 
KOI transits and thus could potentially have been 
found in earlier TPS searches. Regardless, we ex- 
amined all 33 signals more closely. 

For each signal, we phase-folded the light curve 
according to its period and epoch. Close visual 
examination of each signal revealed no credible 
transit-like signal. After checking the original light 
curves, most were determined to be caused by edge 
effects from poorly corrected systematics in the orig- 
inal data. The other signals are spurious for un- 
determined reasons, but could potentially be at- 
tributed to improperly detrended light curves, poor 
transit subtraction, TTVs, or statistical noise. 


5. VISUAL SEARCH FOR PLANETS 
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Figure 3. Phase-folded light curve for all three CPs in Kepler-253. Top: After detrending, but before transit subtraction. 
Bottom: After detrending and transit subtraction. Data points are transparent to emphasize the transit shape. 


Exoplanet systems with multiple transiting plan- 
ets are likelier than other systems to host more dis- 
tant planets that also transit. Some of these will 
only transit once or twice in the Kepler data, which 
are frequently missed by automated search algo- 
rithms. While large planets with 1-2 transits can 
be easy to spot visually, these transits could eas- 
ily be missed or overlooked in the forest of transits 
from known planets. Light curves that are flattened, 
normalized, and have known transits subtracted out 
should then make these planets with 1-2 transits 
stand out more clearly. Therefore, we visually in- 
spected all the transit-subtracted light curves for 
additional transit signals. Three new potential tran- 
sits were found around two stars. 


5.1. Kepler-150 


Two highly significant transits were found in the 
light curve of Kepler-150 and belong to a new 
planet, Kepler-150 f. The transits are visually ap- 


parent in both the PDCSAP flux and the raw SAP 
flux. The second transit slightly overlaps with the 
transit of another planet in the system, but this was 
corrected for in the earlier transit subtraction. A vi- 
sual check confirms that there is no shorter period 
possible in the data. 

The transits were fit with the IDL program TAP 
(Gazak et al. 2012), which is a Markov Chain Monte 
Carlo (MCMC) transit fitter that uses EXOFAST 
(Eastman et al. 2013) to calculate transit models 
(Mandel & Agol 2002) using a wavelet-based likeli- 
hood function (Carter & Winn 2009). TAP fits for 
the basic transit parameters such as the ratio of 
planet radius to stellar radius Rp/R,, the transit 
duration T, the impact parameter b, the midtran- 
sit times, and quadratic and linear limb darkening, 
in addition to white and red noise and a quadratic 
function to correct for improper normalization. A 
circular orbit is assumed. Ten MCMC chains of 
length 200,000 were used to fit the transits. The 
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Figure 4. Phase-folded, flux-normalized light curve for 


Kepler-150 f. Transits by other planets in this window were 
modeled and subtracted out. Blue circles show the first tran- 
sit, while red squares show the second. The black line is the 
best fit from TAP. 


length was picked so as to satisfy the Gelman-Rubin 
statistic (Gelman & Rubin 1992) that tests for non- 
convergence. The planet radius and semi-major axis 
were calculated using the Rp/R, and derived a/R, 
best-fit values from TAP and the stellar radius from 
the KOI catalog cumulative list accessed January 
5, 2017 (Rx = 0.931495 Ro), assuming pseudo- 
Gaussian errors. The reported best-fit values in Ta- 
ble 1 are the median values plus or minus lo error 
bars. The phase-folded, fitted light curve is shown 
in Figure 4. 

Three arguments are used to confirm Kepler-150 f. 
First, the fact that Kepler-150 f is found in a system 
with four other confirmed planets argues strongly 
that it is not a FP. According to Lissauer et al. 
(2012), “almost all of Kepler’s multiple-planet can- 
didates are planets”. They calculate that there 
should be < 1 FP in all systems with 2+ PCs. The 
expected number of FPs for systems with 4+ PCs 
like Kepler-150 would be even lower. 

Secondly, all planets in the Kepler-150 system 
were tested to see if their periods and transit du- 
rations were consistent with orbiting the same star. 
This transit duration analysis has been used pre- 
viously as an additional level of vetting (Steffen 
et al. 2010; Lissauer et al. 2012; Chaplin et al. 2013; 
Cabrera et al. 2014). In a perfectly coplanar, cir- 
cular, edge-on system, and assuming that the plan- 
ets’ radii and masses are much smaller than that 
of the star’s, the transit duration T and the period 
P of the planets are, for any pair of planets, re- 
lated according to the formula T; /Pil _ T;/ iva ' 
where the 7 and 7 indices refer to any two planets 
in the system. The values for Th ca, and Py cae 
were taken from the KOI catalog, while T and Py 
were determined by TAP. The T/P'/? values are 
highly consistent with each other with no planet be- 
ing > 1.6% different from the average value. This 
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Figure 5. Potential single transit in the light curve of 


Kepler-208 at 786.7641 BKJD. Due to its likely nature as a 
pees eclipsing binary FP, we did not attempt a transit 
strongly indicates that all five planets orbit the same 
star and also suggests that they have nearly circular 
orbits. 

Lastly, we performed an analysis with the Python 
program vespa (Morton 2012, 2015; Morton et al. 
2016), which calculates a false positive probabil- 
ity (FPP) for three scenarios of FPs: an eclips- 
ing binary, a hierarchical eclipsing binary, and a 
background eclipsing binary. The vespa analysis of 
Kepler-150 f determines its false positive probability 
to be 0.69%. This does not taken into account ei- 
ther the planet multiplicity argument or the transit 
duration arguments above. 

Combining all three arguments together results in 
a FPP that is < 1%, confirming Kepler-150 f as a 
true exoplanet. It is approximately the size of Nep- 


tune with a planet radius Rp = 3.861929 Reand a 


period P = 637.209470-p185 days. 


5.2. Kepler-208 


Two potential single transits were discovered in 
Kepler-208, a system with four confirmed planets. 
However, one transit has been previously discov- 
ered at Barycentric Kepler Julian Date? BKJD = 
786.7641 by Uehara et al. (2016), who performed 
their own visual checks of KOI systems. The other 
transit is highly suspect. The potential new transit 
is sharply V-shaped (see Figure 5) and has a short 
duration of just JT’ = 2.65 hours. Its morphology 
is most consistent with a single transit of a back- 
ground eclipsing binary, although there is a slight 
possibility that this could be a large, distant planet 
in a glancing transit. However, because it is likely a 
FP, we performed no further analysis of this transit. 


6. DISCUSSION 


3 To convert to Barycentric Julian Date (BJD), use the 
formula BJD = BKJD + 2, 454, 833.0 
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Table 1 
Kepler-150 f properties. 


Parameter 


Best-fit value Unit 


Period (P) 


+0.0164 
637.2047 5.5152 days 


Impact parameter (b) o001s 7% 

a tis : 40.1 
Inclination (7) 89.97 0.3 deg 
Duration (T) res be hours 
Planet radius to stellar radius ratio (Rp/R.) 0.03587) 029 
Semi-major axis to stellar radius ratio (a/R) Sonia 


First midtransit time 


+0.0125 
509.03347* 017, BKJD 


Second midtransit time fidg.o4ne te yee BKJD 
Planet radius (Rp) S86) te Rg 
Semi-major axis (a) LOOe eae AU 
Note. — Best-fit results from the TAP transit fit. Reported values are 


the median and the upper and lower lo error bars. The planet radius and 
semi-major axis were calculated using the stellar radius from the KOI catalog 


(Ry = 0.93170 327 Ro), assuming pseudo-Gaussian errors. Eccentricity was 


held fixed at zero. 


It was not unexpected that these missing planets 
would be rare. There are two reasons why masking 
out data may result in losing planets. One reason is 
that it could reduce the number of detectable tran- 
sits from 3+ to < 3. However, the parameter space 
for which this would occur is small to begin with. 
Planets with periods around 400-500 days (for stars 
with a full ~ 1470 day baseline) are the most sus- 
ceptible to be lost because these are likely to transit 
exactly three times. Due to geometry, planets with 
P = 400 — 500 days are not likely to transit in the 
first place. If such a planet were to transit though, 
then even a single overlapping transit could remove 
it from detectability. However, even an overlapping 
transit would be no guarantee. Long-period plan- 
ets typically have longer durations. Therefore, if 
they overlap with a short-period planet, which is 
the most likely overlapping scenario, then only a 
small portion of the long-period planet’s transit is 
removed, which would allow it to still be detectable. 
In order to remove enough of the long-period tran- 
sit to render it undetectable, either a) it would need 
to overlap with another fairly long-period planet, 
b) two or more short-period planets would need 
to overlap the same long-period transit in different 
spots, or c) a short-period planet and a data gap 
would need to overlap the same long-period transit 
in different spots. Each of these requires an unlikely 
confluence of events. 

A minor confounding factor would be the relative 
impact parameters of the planets with overlapping 
transits. Since planetary systems are usually flat 
with small scatter, inner planets are likelier to have 
lower impact parameters than outer planets in the 
system. Higher impact parameters result in shorter 
durations. Therefore, it is possible, although un- 
likely, that the transit duration of the longer period 
planet would be shorter than (or, at least, compa- 
rable to) that of a shorter period planet in the same 
system. This would make overlapping a larger por- 
tion of the long-period planet’s transit easier and 
thus could more easily render the long-period planet 
undetectable. 


A second way that masking out data may result 
in losing planets is for weak transit signals that 
were on the verge of detectability in the first place. 
The Multiple Event Statistic (MES) is the signal-to- 
noise of a transit signal in the Kepler pipeline (Jenk- 
ins et al. 2002). TPS requires a minimum MES of 
7.1 for detection and classification as a TCE (Jenk- 
ins et al. 2010b; Tenenbaum et al. 2012). Remov- 
ing a full transit or even small portions of a transit 
could result in the MES dropping below this thresh- 
old, thus rendering the planet lost. Subtracting out 
previously found transits rather than removing the 
data points altogether may keep the MES > 7.1. 
On one hand, this might be hard to promote to 
PC or CP status anyway since the FP population is 
dominated by TCEs with three transits and a low 
MES (Mullally et al. 2015). On the other hand, the 
fact that these systems already have 3+ CPs imply 
that these three transit, low MES cases could be 
more easily proven to be planets through statisti- 
cal validation (Lissauer et al. 2012). Note that our 
simulation of planets in Section 2 did not test for 
this. 

A visual search of the data, however, resulted in 
the discovery of Kepler-150 f. This makes Kepler- 
150 just one of 25 stars to host at least five exo- 
planets, according to the NASA Exoplanet Archive 
(Akeson et al. 2013). The demise of the main Ke- 
pler mission, however, has made additional follow- 
up study of this system difficult. A radial veloc- 
ity measurement, for example, would be challeng- 
ing with its long period (P = 637.2 days) and an 
expected radial velocity semi-amplitude of just 1.37 
m/s, assuming a Neptune-mass planet. 


7. CONCLUSION 


A simulation of millions of planets around 114 
stars with 3+ confirmed planets showed that there 
is a low, but non-negligible probability of transiting 
planets being lost after transit masking of known 
KOIs (about 3.3% of transiting planets in the pe- 
riod range of P = 400 — 500 days). We searched 
these same stars for new planetary transits, but in- 
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stead of masking out transits of known KOIs, we fit 
and subtracted them out. However, our search dis- 
covered no credible new transit signals, which was 
consistent with our simulations. 

The original purpose of masking out known in- 
transit data points rather than subtracting them 
out was due to the fact that the subtraction pro- 
cess produced a large number of FPs due to im- 
proper subtraction. This made it impractical to 
do in a time-intensive analysis such as the Kepler 
pipeline despite the risk that it would cause planets 
to be missed. Our analysis, however, demonstrates 
the validity and robustness of the Kepler pipeline’s 
choice to use transit masking over transit subtrac- 
tion. 

However, a visual follow-up of the transit- 
subtracted light curve revealed the existence of 
a Neptune-sized exoplanet, Kepler-150 f (Rp = 
3.86 Re), making Kepler-150 one of the few stars 
with 5+ known planets. Because of its long period 
(P = 637.2 days), only two transits are contained in 
the data, which made it undetectable to the Kepler 
pipeline. The authors contribute its discovery to 
the subtraction of known planet transits from the 
light curve. This discovery suggests the possibil- 
ity that improved light curve flattening and transit 
subtraction, or simply better eyes, may result in the 
discovery of new, long period exoplanets. 
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Table 2 
Stars with 3+ confirmed used in this study. 


Kepler KIC Number of confirmed Number of KOI Name of non-fitted 
name planets (not fit) candidates (not fit) planets and KOIs 
Kepler-11 6541920 6 0 
Kepler-18 8644288 3 0 
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Table 2 — Continued 


Kepler KIC Number of confirmed Number of KOI Name of non-fitted 
name planets (not fit) candidates (not fit) planets and KOIs 

Kepler-20 6850504 5 0 
Kepler-30 3832474 3 0 
Kepler-31 9347899 3 Al 
Kepler-33 9458613 5 0 
Kepler-37 8478994 3 0 
Kepler-42 8561063 3 0 
Kepler-48 5735762 3 0 
Kepler-49 5364071 4 0 
Kepler-52 11754553 3 0 
Kepler-53 5358241 3 0 
Kepler-54 7455287 3 0 
Kepler-55 8150320 5 (1) 0 Kepler-55 c 
Kepler-58 4077526 3 1 
Kepler-60 6768394 3 0 
Kepler-62 9002278 5 0 
Kepler-65 5866724 3 0 
Kepler-79 8394721 4 0 
Kepler-80 4852528 4 a 
Kepler-81 7287995 3 0 
Kepler-82 7366258 4 0 
Kepler-83 7870390 3 0 
Kepler-84 5301750 5 0 
Kepler-85 8950568 4 0 
Kepler-89 6462863 4 0 
Kepler-102 10187017 5 0 
Kepler-104 6678383 3 0 
Kepler-107 10875245 4 0 
Kepler-114 10925104 3 0 
Kepler-122 4833421 5 0 
Kepler-124 11288051 3 0 
Kepler-127 9451706 3 0 
Kepler-130 5088536 3 0 
Kepler-132 6021275 3 1 
Kepler-138 7603200 3 0 
Kepler-142 10982872 3 0 
Kepler-149 3217264 3 0 
Kepler-150 5351250 4 0 
Kepler-164 10460984 3 1 (1) KOI 474.03 
Kepler-169 5689351 5 0 
Kepler-171 6381846 3 0 
Kepler-172 6422155 4 0 
Kepler-174 8017703 3 0 
Kepler-176 8037145 3 1 
Kepler-178 9941859 3 0 
Kepler-184 7445445 3 0 
Kepler-186 8120608 5 0 
Kepler-194 10600261 3 0 
Kepler-197 12068975 4 0 
Kepler-203 6062088 3 0 
Kepler-206 6442340 3 0 
Kepler-207 6685609 3 0 
Kepler-208 7040629 4 0 
Kepler-215 8962094 4 0 
Kepler-219 9884104 3 0 
Kepler-220 9950612 4 0 
Kepler-221 9963524 4 0 
Kepler-222 10002866 3 0 
Kepler-223 10227020 4 (1) 0 Kepler-223 e 
Kepler-224 10271806 4 0 
Kepler-226 10601284 3 0 
Kepler-228 10872983 3 0 
Kepler-229 10910878 3 0 
Kepler-235 4139816 4 0 
Kepler-238 5436502 5 0 
Kepler-244 6849310 3 0 
Kepler-245 6948054 3 (2) 1 Kepler-245 b, Kepler-245 d 
Kepler-249 7907423 3 0 
Kepler-250 8226994 3 0 
Kepler-251 8247638 4 0 
Kepler-253 8689373 3 0 
Kepler-254 9334289 3 0 
Kepler-256 9466668 4 0 
Kepler-257 9480189 3 0 
Kepler-265 5956342 4 0 
Kepler-267 10166274 3 0 
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Table 2 — Continued 


Kepler KIC Number of confirmed Number of KOI Name of non-fitted 
name planets (not fit) candidates (not fit) planets and KOIs 


Kepler-272 10426656 3 
Kepler-275 3447722 3 
Kepler-276 3962243 3 
Kepler-282 8609450 4 
Kepler-286 10858691 4 
Kepler-288 4455231 3 
Kepler-292 6962977 5 
Kepler-295 9006449 3 (1) 
Kepler-296 11497958 5 (1) 
Kepler-298 11176127 3 
Kepler-299 11014932 4 
Kepler-301 11389771 3 
Kepler-304 5371776 3 
Kepler-305 5219234 3 
Kepler-306 5438099 4 
Kepler-310 10004738 3 
Kepler-325 9471268 3 
Kepler-327 8167996 3 
Kepler-331 4263293 3 
3 
3 
3 
4 
3 
4 
( 
3 
3 
( 
3 
3 
3 
3 
4 
5 
3 
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Be 


SPOOF ONTDOOCOCOORFCAOOCOCOCOCOOCOCOOFRFROCOCOCCOCOCOCCCO7O 


KOI 1198.04 


= 
Vm 


Kepler-295 d 
Kepler-296 e 


Kepler-332 10328393 
Kepler-334 10130039 
Kepler-336 6037581 
Kepler-338 5511081 
Kepler-339 10978763 
Kepler-341 7747425 
Kepler-342 9892816 3 
Kepler-350 4636578 
Kepler-354 6026438 
Kepler-357 8164257 
Kepler-363 6021193 
Kepler-372 11401767 
Kepler-374 6871071 
Kepler-399 5480640 
Kepler-402 7673192 
Kepler-444 6278762 
Kepler-445 9730163 
Kepler-446 8733898 


= 
VY 


w 
= 
ro, 


Kepler-357 d 


Note. — Stars included in this study, each with 3+ CPs. 
Some also have additional PCs. The column “Number of 
confirmed planets (not fit)” refers to the number of CPs in 
that system, while the number in the parentheses, if appli- 
cable, are how many of those CPs were unable to be fit. 
The column “Number of KOI candidates (not fit)” is simi- 
lar, but for PCs that have not been confirmed. The names 
of the non-fitted KOIs are in the last column. 


